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A  SIMPLE  MODEL  OF  AEROSOL  OPTICAL  PROPERTIES  IN  A  WELL-MIXED 
MARINE  BOUNDARY  LAYER  FOR  A  WAVELENGTH  OF  0.50^ 

1  INTRODUCTION  * 

It  is  well  known  that  absorption  and  scattering  by 
aerosol  particles  in  the  atmosphere  has  an  important  effect 
on  radiative  transfer  and  the  radiation  budget  of  the  earth. 
The  aerosol  optical  thickness,  being  the  vertically 
integrated  value  of  extinction  due  to  aerosols,  is  thus  an 
important  parameter  affecting  climate.  Knowledge  of  the 
vertical  profile  of  aerosol  extinction  is  also  of  importance 
in  studies  concerned  with  the  performance  of  the  numerous 
electrooptical  systems  now  in  use.  These  include 
communications  and  remote  sensing  systems  and  a  variety  of 
systems  having  military  applications. 

Since  the  oceans  cover  roughly  three- fourths  of  the 
earth's  surface,  knowledge  of  aerosol  extinction  in  the 
marine  boundary  layer  is  of  obvious  importance.  An 
inversion  capped  convective  boundary  layer  is  a  fairly 
common  and  often  persistent  feature  of  the  lower  atmosphere 
over  many  areas  of  the  ocean.  The  inversion  marks  the 
boundary  between  a  turbulent  well -mixed  layer  below  and  a 
stably-stratified,  non- turbulent  free  atmosphere  above.  The 
structure  and  evolution  of  the  marine  boundary  layer  is 

usually  modeled  by  means  of  mixed  layer  models  (Stage  and 
Blaniueript  approvad  Fabruary  20, 1986. 


Businger,  1981;  Tennekes  and  Driedonks,  1981;  Wakefield  and 
Schubert,  1981).  According  to  these  models,  the  region 
between  the  top  of  the  surface  layer  (~10  m)  and  the  base  of 
the  inversion  is  assumed  to  be  fully  mixed  with  constant 
values  of  quantities  that  are  conserved  during  adiabatic 
expansion.  Davidson  et  al.  (1984)  relaxed  the  zero 
gradient  assumption  to  permit  a  small  constant  gradient  of 
conservative  quantities.  For  a  dry  convective  boundary 
layer  the  conserved  properties  are  potential  temperature, 
water  vapor  mixing  ratio  and  aerosol  number  concentration 
(particles  per  gram  of  air),  while  for  a  cloud-topped  mixed 
layer  the  conserved  quantities  are  the  equivalent  potential 
temperature  and  the  total  water  mixing  ratio. 

This  paper  describes  a  model  of  the  aerosol  optical 
thickness  and  the  vertical  profile  of  aerosol  extinction  for 
the  case  of  a  well -mixed  marine  boundary  layer.  Presented 
in  the  form  of  both  nomograms  and  analytical  approximation 
formulas,  the  model  predicts  the  vertical  variation  of  the 
aerosol  extinction  coefficient  in  the  boundary  layer  given 
the  relative  humidity,  temperature  and  extinction 
coefficient  at  the  surface  (shipboard  level)  and  the  height 
of  the  boundary  layer.  Conversely,  given  the  relative 
humidity  and  extinction  coefficient  at  the  surface  and  an 
estimate  of  the  optical  thickness  of  the  boundary  layer,  the 
model  gives  an  estimate  of  the  depth  of  the  mixed  layer. 
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Although  the  model  is  designed  for  the  cloud* free  boundary 
layer  case,  it  should  be  quite  applicable  to  a  boundary 
layer  containing  scattered  fair  weather  cumulus  as  long  as 
the  fractional  cloud  cover  is  less  than  a  few  tenths. 


3 


where  n(r)  is  the  number  of  particles  per  cm  per  micron 
radius  interval  and  is  the  extinction  efficiency  factor 
computed  from  Mie  theory  for  spherical  homogeneous 
particles.  Since  we  are  assuming  well  mixed  conditions,  the 
change  in  n(r}  with  altitude  is  due  mainly  to  an  increase  in 
relative  humidity  with  only  a  small  effect  due  to  the 
decrease  in  particle  concentrations  as  a  result  of  the 
decrease  in  air  density  with  height. 

Vertical  profiles  of  relative  humidity  for  a  well -mixed 
layer  are  shown  in  Fig.  1  for  various  values  of  surface 
(i.e.,  10  ffl)  relative  humidity.  A  surface  temperature  of 
75°F  was  assumed.  Relative  humidity  is  defined  as  the 
quotient  of  the  observed  water  vapor  mixing  ratio  to  the 
saturation  mixing  ratio: 


Percent  relative  humidity  (RH)  is  S  times  100.  The 
saturation  mixing  ratio  is  a  function  of  both  temperature 
and  pressure  and  is  defined  as 


m 


s 


0.622es 
P  *  *3  ’ 


where  e,  is  the  sacuration  vapor  pressure  and  P  is 


atmospheric  pressure. 


In  computing  the  relative  humidity 


profiles  shown  in  Fig.  1  it  was  assumed  that  m  is  constant 
with  height,  that  temperature  decreases  at  the  dry  adiabatic 
lapse  rate,  -  -9.76°C/km  and,  therefore,  that  pressure 
changes  according  to 


where  and  are  the  pressure  and  temperature  at  the 
surface  and  a  -  Here  g  is  the  acceleration 
of  gravity  and  is  the  gas  constant  of  dry  air. 

Vertical  profiles  of  the  normalized  extinction 
coefficient,  ,  were  computed  from  the  relative 
humidity  profiles  in  Fig.  1  and  a  model  of  aerosol 
extinction  as  a  function  of  relative  humidity.  The 
dependence  of  aerosol  extinction  on  relative  humidity  was 
calculated  from  Eq.  (3)  using  the  model  of  particle  size 
versus  relative  humidity  shown  in  Fig.  2  and  the  marine 
aerosol  size  distribution  shown  in  Fig.  3.  This  size 
distribution  is  a  smooth  fit  to  a  size  distribution  that  was 
measured  approximately  600  km  off  the  east  coast  of  the  U.S. 
on  March  13,  1983  during  a  cruise  of  the  USNS  Lynch. 
Particles  smaller  than  O.Sfim  were  measured  with  the  NRL 
differential  mobility  analyzer  while  particles  larger  than 
0.6(xm  were  measured  with  a  gelatin  replication  impactor 
(Calspan  Corp.).  The  relative  humidity  at  the  time  of  the 
measurements  was  52%. 
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Fig.  2  shows  the  ratio  of  particle  radius  at  a  given 
humidity  to  particle  radius  in  dry  state  plotted  as  a 
function  of  relative  humidity.  Particles  smaller  than  0.5|j.m 
are  assumed  to  grow  according  to  curve  1  while  particles 
larger  than  1.5|j.m  are  assumed  to  follow  growth  curve  2.  The 
growth  curve  of  intermediate  sized  particles  is  assumed  to 
lie  between  these  two  curves.  Growth  curve  1  assumes  that 
particles  are  90%  soluble,  the  soluble  material  being  a 
mixture  of  ammonium  sulfate  and  ammonium  nitrate.  It  is 
based  on  the  measurements  of  Hanel  and  Lehmann  (1981)  of  the 
water  uptake  as  a  function  of  relative  humidity  by  samples 
of  atmospheric  aerosol  material  that  contained  these  two 
salts  as  the  primary  water-soluble  constituents.  This 
growth  model  for  the  smaller  marine  aerosol  particles  is 
consistent  with  recent  measurements  (Hoppel  et  al.,  1985)  of 
the  volatility  of  submicron  particles  over  the  Pacific.  In 
these  experiments,  particles  were  subjected  to  a  level  of 
heating  that  would  evaporate  the  more  volatile  aerosol 
components  without  affecting  the  highly  involatile  sea-salt 
component.  It  was  found  that  nearly  all  the  particles  were 
highly  volatile,  indicating  that  they  likely  were  composed 
of  gas-phase  conversion  products  such  as  sulfates  and 
nitrates.  Curve  2,  describing  the  growth  of  the  larger 
particles,  is  the  mean  of  Hanoi's  (1976)  sea  spray  aerosol 
growth  model  (Hanoi's  model  2)  and  his  model  3  for  maritime 


averaged  the  growth  curves  for  increasing  and  decreasing 
relative  humidity. 

In  computing  the  growth  curves,  the  effect  of  particle 
curvature  on  the  equilibrium  relative  humidity  was  neglected 
so  that  the  curves  are  strictly  valid  for  infinitely  large 
particles.  The  error  in  r/r^  due  to  the  neglect  of  particle 
curvature  is  less  than  about  5%  for  r^  >0.lM.m  and  RH<97%. 
Furthermore,  since  particles  smaller  than  O.lfxm  contribute 
less  than  a  few  percent  to  the  extinction  coefficient  at 
Q.5|xm  for  typical  marine  aerosol  size  distributions,  the 
resulting  error  in  the  computed  variation  of  aerosol 
extinction  with  relative  humidity  is  quite  small. 


Fig.  4  shows  the  calculated  variation  of  aerosol 
extinction  with  relative  humidity  for  a  wavelength  of 
O.SO^m.  The  results  are  presented  as  a  plot  of  the 
normalized  extinction  coefficient  as  a  function  of  relative 
humidity.  The  dependence  of  the  normalized  extinction  on 
relative  humidity  is  determined  primarily  by  the  aerosol 
growth  model  assumed  (i.e.,  on  aerosol  composition)  and 
depends  only  weakly  on  the  particle  size  distribution.  We 
computed  the  normalized  extinction  for  a  number  of  size 
distributions  measured  on  the  1983  cruise,  including  a 
typical  distribution  for  the  remote  tropical  Atlantic,  and 
found  a  maximum  variation  of  p^(RH}/p^(RH>52%)  of  about  20% 
at  97%  relative  humidity  with  a  smaller  variation  at  lower 
humidities . 


From  the  variation  of  relative  humidity  with  height 
(Fig.  1)  and  the  variation  of  extinction  with  relative 
humidity  (Fig.  4)  we  computed  the  profiles  of  normalized 
aerosol  extinction  shown  in  Fig.  5.  The  abscissa  is  the 
aerosol  extinction  at  height  z  divided  by  the  aerosol 
extinction  at  the  surface  (10-m  level).  Each  curve  is  for  a 
different  value  of  surface  relative  humidity  as  identified 
in  the  key.  Thus,  if  we  know  the  aerosol  extinction 
coefficient  and  relative  humidity  at  the  surface  and  the 
height  of  the  boundary  layer,  then  Fig.  5  can  be  used  to 
obtain  the  profile  of  aerosol  extinction  within  the  boundary 
layer.  The  accuracy  of  the  predicted  profile  of  aerosol 
extinction  depends  on  the  accuracy  to  which  the  extinction 
coefficient  and  relative  humditity  at  the  surface  can  be 
determined.  The  aerosol  extinction  coefficient  at  the 
surface  can  be  obtained  from  direct  measurements  of  the 
scattering  or  extinction  coefficient,  from  visibility 
measurements  via  the  Roschmieder  formula  (i.e.,  V-3.91/Pg) 
or  from  Nie  calculations  using  a  predicted  or  measured  sur¬ 
face  aerosol  size  distribution. 

The  profiles  of  the  normalized  aerosol  extinction 
coefficient  in  Fig.  5  were  computed  only  up  to  the  height  at 
which  the  relative  humidity  reached  97%.  There  are  two 


principal  reasons  why  the  profiles  are  truncated  at  this 


height.  First,  since  the  vertical  resolution  of  temperature 
and  humidity  profiles  obtainable  with  standard  radiosondes 
is  only  about  150  feet,  the  height  of  the  boundary  layer  can 
only  be  resolved  to  within  150  feet  or  so.  Within  a 
distance  of  150  feet  the  relative  humidity  would  increase 
from  97%  to  almost  100%  and  the  extinction  would  change  by  a 
factor  of  about  three.  Thus,  in  cases  when  the  relative 
humidity  at  the  top  of  the  cloud-free  (or  mostly  cloud  free) 
boundary  layer  is  above  97%,  uncertainties  in  determining 
the  height  of  the  boundary  layer  can  lead  to  large  errors  in 
the  predicted  aerosol  extinction  coefficient.  Second,  if 
the  boundary  layer  is  deep  enough  for  the  relative  humidity 
to  exceed  97%  then  there  is  a  good  chance  that  the  boundary 
layer  will  be  cloud  covered  and  that  the  optical  properties 
in  the  region  between  the  97%  humidity  level  and  cloud  base 
will  be  affected  by  the  presence  of  the  clouds  (e.g., 
through  turbulent  mixing  at  cloud  base). 

Eq.  (2)  shows  that  the  aerosol  optical  thickness  of  the 
boundary  layer  can  be  expressed  as  the  product  of  the 
extinction  coefficient  at  the  surface  and  the  vertically 
integrated  value  of  Pg( z )/Pg( z^) .  The  integral  in  Eq.  (2) 
was  evaluated  using  the  profiles  of  aerosol  extinction  given 
in  Fig.  5.  Fig.  6  is  a  plot  of  this  integral  as  a  function 
of  boundary  layer  height,  H,  for  various  values  of  surface 
relative  humidity.  The  key  in  Fig.  5  identifies  the  values 


»**  V'"  ^ 
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of  surface  relative  humidity  corresponding  to  each  curve. 
The  dashed  curve  in  Fig.  6  gives  the  value  of  the  integral 
if  there  were  no  change  in  aerosol  extinction  with  height. 
In  this  case  the  integral  simply  has  the  value  H  —  .  To 
estimate  the  optical  thickness  of  the  boundary  layer  we  use 
Fig.  6  to  obtain  the  appropriate  value  of  the  integral  in 
Eq.  (2)  and  then  multiply  this  value  by  the  aerosol 
extinction  coefficient  at  the  surface. 


3  ANALYTICAL  MODEL 

Although  the  nomograms  (Figs.  5  and  6)  are  very  useful, 
interpolation  is  required  for  other  values  of  surface 
relative  humidity.  It  is  clearly  advantageous  to  have  an 
analytical  form  of  the  model  to  facilitate  computer  cal¬ 
culations  of  the  aerosol  extinction  profile.  One  such  model 
is  developed  below. 

Relative  humidity  may  be  calculated  from  temperature  T 
and  dew  point  T^  by  means  of  the  expression 


"l9.5  (Tj  -  T) 
269.9  +  T  +  Tj 


where  S  is  the  saturation  ratio  and  T  and  T^  are  in  ®C. 
This  expression  is  identical  to  the  one  developed  by  Shettle 
(1978)  except  that  in  Shettle's  original  expression  the 


coefficient  of  T, -T  in  the  numerator  is  19.772.  Shettle's 

a 

approximation  has  a  RMS  error  of  0.11%  for  0 < (T-T^ )  < 40°C. 
Near  the  ocean  surface  one  is  normally  dealing  with  a 
smaller  range  of  T*T^ . 

Denoting  the  values  of  T  and  at  shipboard  height 

( z»Zq )  by  Tg  and  T^g  we  have,  assuming  adiabatic  mixing. 


T(Z)  -  Tg  +  (Z  -  Zg) 


^d'^*  ■  ■'do  'd.'^  -  "d>- 


where  and  F^^  are  the  adiabatic  lapse  rates  of  tempera¬ 
ture  and  dew  point,  respectively.  We  take  F^^  »  -.00175°C/m 
and  Fa  -  -.  00976 ®C/m. 

Combining  Eqs.  (6)  and  (7)  yields 


S(z)  ■  exp 


19.5  (T.  -  T  )  ■*■  .156  (r-z  ) 

do  o  o 

269.9  +  T  +  T .  -  .0115  (z-z  ) 

o  do  o 


where  S(z}  is  the  saturation  ratio  at  height  z.  Eq.  (8) 
gives  the  relative  humidity  at  height  z  as  a  function  of  the 
temperature  and  dew  point  at  the  surface.  The  surface  dew 
point  temperature  can  be  computed  from  the  relative  humidity 
and  temperature  at  the  surface  using 


19. 5T  +  (269. 9+T  )  InS 


/Q'i 


where  S  is  the  saturation  ratio  at  the  surface. 


The  dependence  of  aerosol  extinction  on  relative 
humidity  given  in  Fig.  4  can  be  approximated  by  the 
expression 

8^(S)  r  S  -  .52 

_S -  .  I  + - 

8^(S-.52)  L  1-017  -  S 


This  equation  is  accurate  to  within  10%  in  the 
range  0.73  <3  4:0.82  and  to  within  1%  to  5%  for  values  of  S 
outside  this  range.  The  ratio  of  extinction  at  relative 
htimidity  Sj  to  that  at  relative  humidity  ,  where 


S2>Sj^^0.52,  may  be  written 


VV  ,  °.«2>  / 

9,(Sj)  •  a^<S..32)/  3^(S-.52) 


which  reduces  to 


w 

w 


(11) 


Substituting  S(z)  for  $2  and  Sg  for  in  Eq.  (11) 


yields 


Profiles  of  aerosol  extinction  in  the  marine  boundary  layer 
can  be  computed  from  Eqs.  (8),  (9)  and  (12)  provided  we  know 
the  height  of  the  boundary  layer,  and  the  temperature, 
relative  humidity  and  extinction  coefficient  at  the  surface. 

The  extinction  coefficient  at  height  z  can  be  adjusted 
for  the  decrease  in  particle  number  density  (particles  per 
cm^)  with  height  resulting  from  the  decrease  in  air  density 
by  multiplying  £he  right  side  of  Eq.  (12)  by  p(z)/p(Zo)r 
which  is  the  ratio  of  the  density  of  air  at  height  z  to  the 
density  of  air  at  the  surface.  The  change  in  air  density 
with  altitude  may  be  approximated  by 

p(z)/P(Zo)  ■ 

where  A  -  7.94  x  10"^  m'l. 

4  synnARY 

A  simple  model  to  predict  the  vertical  profile  of 
aerosol  extinction  and  the  aerosol  optical  thickness  in  the 
marine  boundary  layer  is  presented.  The  model  assumes  that 
the  boundary  layer  is  well  mixed,  so  that  the  potential 
temperature,  water  vapor  mixing  ratio  and  dry  aerosol  size 
distribution  are  constant  with  height  from  the  10 -m  level 
(shipboard  level)  to  the  inversion.  Required  inputs  to  the 
model  are  the  temperature,  relative  humidity  (or  dew  point) 
and  extinction  coefficient  at  the  surface  and  the  height  of 


the  boundary  layer.  Particles  larger  than  1.5(jim  radius  are 
assumed  to  grow  with  relative  humidity  according  to  curve  2 
in  Fig.  2.  This  curve  is  a  mean  between  Hanel's  (1976) 
growth  model  2  (sea  spray  aerosols)  and  model  3  (maritime 
aerosol  over  the  Atlantic).  Particles  smaller  than  O.Spjn 
follow  growth  curve  1  which  describes  the  behavior  of 
particles  which  are  90%  soluble,  the  soluble  component 
consisting  of  ammonium  sulfate  and  ammonium  nitrate. 
Intermediate  sized  particles  have  a  growth  curve  between 
curves  1  and  2,  reflecting  a  transition  from  particles 
composed  primarily  of  gas*to*particle  conversion  material  to 
sea  salt  particles. 

The  model  is  presented  in  two  forms:  (1)  Nomograms 
(Figs.  5  and  6)  which  give  the  ratio  of  the  aerosol 
extinction  coefficient  at  height  z  to  aerosol  extinction  at 
the  surface  (Fig.  5)  and  the  vertically  integrated  value  of 
aerosol  extinction  as  a  function  of  boundary  layer  height 
and  surface  relative  humidity  (Fig.  6)  and  (2)  Analytical 
approximation  formulas  [Eqs.  (8),  (9)  and  (12)]  from  which 
one  can  compute  aerosol  extinction  as  a  function  of  height 
given  the  required  input  values. 
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RELATIVE  HUMIDITY  (%) 

Models  of  particle  size  as  a  function  of  relative 
humidity  used  in  the  calculations.  Curve  1  is 
for  particles  smaller  than  0.5|im  radius  while 
curve  2  is  for  particles  larger  than  l.Sum. 


with  particle  growth  curves  to  compute  aerosol 
extinction  as  a  function  of  relative  humidity. 
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Ratio  of  the  aerosol  extinction  coefficient  at 
height  z  to  aerosol  extinction  at  the  surface, 
plotted  as  a  function  of  z.  Each  curve  is  for  a 
different  value  of  surface  relative  humidity  as 
identified  in  the  key. 
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The  value  of  the  integral  in  Eq.  (2)  plotted  as  a 
function  of  boundary  layer  height,  H,  for  various 
values  of  surface  relative  humidity.  [See  Fig.  5 
for  a  key  to  the  curves]. 


REFERENCES 


Davidson,  K.  L.,  C.  W.  Fairall,  P.  Jones  Boyle,  and  G.  E. 
Schacher,  1984:  Verification  of  an  atmospheric 

mixed-layer  model  for  a  coastal  region.  J.  Climate 
Aopl .  Meteor . ,  23 ,  617-636. 

Hanel,  G.,  1976:  The  properties  of  atmospheric  aerosol  par¬ 
ticles  as  functions  of  the  relative  humidity  at  thermo¬ 
dynamic  equilibrium  with  the  surrounding  moist  air. 
Advances  in  Geophysics,  Vol.  19,  Academic  Press, 

73-188. 

Hanel,  G.,  and  M.  Lehmann,  1981:  Equilibrium  size  of  aerosol 
particles  and  relative  humidity:  New  experimental  data 
from  various  aerosol  types  and  their  treatment  for 
cloud  physics  applications.  Contrib.  Atmos.  Phys.,  5£, 
57-71. 

Hoppel,  W.  A.,  G.M,  Frick,  and  R.E.  Larson,  1985:  Effect  of 
nonprecipitating  clouds  on  the  aerosol  size  dis¬ 
tribution  in  the  marine  boundary  layer.  Geophvs.  Res. 
Letters  (In  Press). 

Shettle,  E.  P.,  1978;  Analytic  expressions  for  the  atmos¬ 

pheric  water  vapor  content.  OPAQUE  N/US-7810,  p.  5. 

Stage,  S.  A.,  and  J.  A.  Businger ,1981 :  A  model  for  entrain¬ 
ment  into  a  cloud- topped  marine  boundary  layer.  Part 
I:  Node!  description  and  application  to  a  cold-air 

outbreak  episode.  J.  Atmos.  Sci..  38 ,  2213-2229. 


Tennekes,  H.,  and  G.  N.  Driedonks,  1981:  Basic  entrainment 
equations  for  the  atmospheric  boundary  layer. 
Bound. -Laver  Meteor..  20 .  515-531. 

Wakefield,  J.S.,  and  W.  H.  Schubert,  1981:  Nixed-layer  model 
simulation  of  eastern  North  Pacific  stratocumulus . 
Mon.  Wea.  Rev..  109,  1952-1968. 


